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Abstract. In this paper we construct examples of commutative rings of difference 
operators with matrix coefficients from representation theory of quantum groups, 
generalizing the results of our previous paper [ES] to the q-deformed case. A general- 
ized Baker- Akhiezer function \? is realized as a matrix character of a Verma module 
and is a common eigenfunction for a commutative ring of difference operators. 

In particular, we obtain the following result in Macdonald theory: at integer 
values of the Macdonald parameter k, there exist difference operators commuting 
with Macdonald operators which are not polynomials of Macdonald operators. This 
result generalizes an analogous result of Chalyh and Veselov for the case q = 1, 
to arbitrary q. As a by-product, we prove a generalized Weyl character formula for 
Macdonald polynomials (= Conjecture 8.2 from [FV]), the duality for the ^-function, 
^{■^ . and the existence of shift operators. 
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Let iV be a positive integer. Let be the algebra over the field C(q) gener- 
ated by the field of rational functions C(q, Xi, Xn) and commuting operators 
T x 1 , ...jTjy 1 , with commutation relations 

Ti o f(q, Xi, ...Xi, X N ) = f(q, Xi, qX h ...,X N ) o Ti. 

This algebra is called the algebra of (/-difference operators in n variables with ra- 
tional coefficients. Elements of this algebra are called difference operators. 

Let V be a finite-dimensional vector space over C. Introduce the algebra Dg (V) 
of difference operators with matrix coefficients 

Dq (V) = <8> End(V) 

Let q be a simple finite-dimensional Lie algebra over C of rank r, and let U q g be 
the corresponding quantum group. In [EK], to any finite dimensional representation 
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U of U q Q was assigned a family of commuting difference operators D c parametrized 
by Weyl group invariant trigonometric polynomials c on the Cartan subalgebra f) 
of q. These operators are constructed as follows. 

Let M\ be the Verma module over U q Q with highest weight A and highest weight 
vector v\. Let U[0] be the zero weight subspace of U. For any u G U[0], define the 
intertwining operator <£>^ : M\ — > M\®U by the condition &™v\ = v\®u+J2 Wi®Ui, 
where Wi are homogeneous vectors of weights \ii < A. This operator is defined 
for generic A. For any weight u, let Proj| M ^^j: M\ — > M\ be the homogeneous 

projector to the subspace M^[z/] of weight v. Let ty\{X\, X r ) be the function 
with values in End(U[0]) such that for any u £ U[0] 

^(X 1 ,...,X r ) W = ^Xr\..^Tr| MA (Proj| MAN o^oProj| MAM ) 

Proposition 1.1. [EK] For any Weyl group invariant function c(A) on h* of the 
form 

(1-1) c(X) = J2^q 2M , c,EC(q) 

A*eP 

there exists a unique difference operator D c e £>g(£/[0]) such that 

D c 4>\ = c(A + p)$x 

For any root a of 0, let k a = ma,x{n\U[na] ^ 0}, where U[n] is the subspace of 
weight fi in U. Let R(U) be the ring of functions on h* of the form (1-1) such that 
for any positive root a of q 

c(^A-— J=c(^A + — J, n=l,...,k a 

whenever (A, a) = 0. The main result of this paper is the following theorem, proved 
in Chapter 4 of this paper. 

Theorem 1.2. There exists an infective homomorphism £ : R(U) — > Dg(?7[0]) 
such that for any Weyl group invariant element c G R one has £(c) = D c . For any 
c G R(U), the operator £(c) is defined by the equation 

£(c)^ A = C(A + p)^ A - 

We will denote £ (c) by _D C for any £ G R. 

In the case when g = sIn (type Ajv-i), we can choose representation [7 to be 
5 ,fc7V l/ 7 where V is the fundamental representation, in which case the space U[0] 
is 1-dimensional. Then the operators D c for symmetric functions c are conjugate 
to Macdonald operators, corresponding to t = q k+1 . Namely, if q are elementary 
symmetric functions, then {D Cl } are simultaneously conjugate to 

lCl,...,N,\I\=li€l,j£I 1 3 iei 

(in suitable coordinates). In this case, the numbers k a are all equal to k; so we will 
denote the algebra R(U) by Rk- From Theorem 1.2 we get (See Chapter 5): 
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Theorem 1.3. For any positive integer k, there exists an infective homomorphism 
£ : Rk — > such that = Mi, I = 1, . . . , N. The function ifj\ is a common 
eigenfunction of the operators ((c),c£ Rk with eigenvalue c(A + p). 

Note that Theorem 1.3 is a special property of Macdonald's operators at integer 
values of k. If k is not an integer, one can show that the centralizer of Aii, ■ ■ ■ , M-n 
in reduces to the polynomial algebra of A4±, . . . , M.n- We call this special prop- 
erty at integer values of k "algebraic integrability of Macdonald operators" , by anal- 
ogy with the case differential operators which was treated in [CV1,CV2,VSC,ES]. 
In this sense, the results of this paper is precisely a (/-deformation of the results of 
[ES]. 

As a by-product, we obtain several results in Macdonald's theory. Namely, we 
prove the partial Weyl group symmetry of the -^-function, a generalized Weyl char- 
acter formula for Macdonald's polynomials (which coincider with Conjecture 8.2 
in [FV]), an explicit formula for the -^-function in terms of shift operators, and 
symmetry of the -^-function with respect to the interchange A «-> x. 

The paper is organized as follows. In Section 2 we recall basic facts about repre- 
sentations of quantum groups and intertwining operators. In Section 3 we introduce 
the ^-function as matrix trace of an intertwining operator, and prove its proper- 
ties. In Section 4 we explain how to construct a commutative ring of difference 
operators from the \l/-function. In Section 5 we review some facts from Macdonald 
theory for root system A n and explain how to obtain them from our construction. 
In Appendix we show how our technique works in the simplest example. 
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2. Quantum groups and their representations 

Notation. Let q be a simple (finite-dimensional) complex Lie algebra of rank r 
with fixed diagonalizable Cartan matrix A = (a^), i, j = 1, . . . , r, and let di, . . . , d r 
be positive relatively prime integers such that the matrix B = (6^) = {dittij) is 
symmetric. We denote its Cartan subalgebra by h. Let ctj e f)*, i = 1, . . . , r denote 
simple roots, R - the corresponding root system, R + and R~ - the sets of positive 
and negative roots, respectively. 

The invariant Killing form (■, •) on {)* is defined by (ctj, aj) = dittij. Let Ai, . . . , A r e 
()* be fundamental weights, i.e. (Aj, aj) = Sij, i,j = 1, . . . , r. Put p = Yli=i ^-i- 
Denote 

q = Q+ = z + a " P = J2 ZA " p + = Yl Z + A - 

For /i, v G P we write n>v if p — v <E Q+. 

Let W be the Weyl group of g. The Weyl group generators Sj act on h* by simple 
root reflections 

Si- p = p- 2- 

{ai, ai) 

We also introduce a shifted action of Weyl group by 



w p ■ p = w(p + p) — p 
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For w G W let l(w) denote the length of w, i.e. the number of generators in a 
reduced decomposition w = ■ . . . ■ Sj r 

Quantum groups. The quantum group U q Q, associated to simple Lie algebra 0, 
is a Hopf algebra over C(q) with generators E t , F t , i = 1, . . . , r and relations: 

KiKj = KjKi, K.F.j = <{',■ F. :i K,. K , /•', = q- " : !■)!<, 



EF — F E — <$■■ 



qi-q, 



-i ' 



l-a,- 



Em: 



k=0 
1 — aa 



fc=0 



1 — a 
fc 

1 — a 
fc 



>3 



■£,1 a,ij k jj, jjik 



F- aij k F 3 F* = 0, i + j. 



where qi = q 1 and we used notation 

! 



n 



n r , — n 



-, [n] q \ = [l] q -[2] q -....[n] q , [n] q = 



q — q 



, [k\ q \[n-k\ q \ q-q- 

Comultiplication A, antipode S and counit e in U q Q are given by 

A(Ei) = E, <g> 1 + K % <g> Ei, A(Fi) = F, <g> K~ l + 1 <g> F t , A(K t ) =K i ®K i 

S(E t ) = —K~ 1 E i , S(F t ) = —FiKi, S(K t ) = K~ x 
e{Ei) = 6(F) = 0, e(Ki) = 1 
We define a C-algebra involution uj of U q Q by 

to(E i ) = -F i , u(Fi) = -Ei, u(Ki) = Ki, u(q) = q' 1 

We have a decomposition of vector spaces U q Q = U~®U°®U + , where U~ (resp. U + ) 
is the subalgebra generated by Fi (resp. Ei), and U° is generated by K i7 K~ x , i = 
l,...,r. 

Verma modules. For any A G f)* we can introduce Verma module M\ over U q Q, 
i.e. W~-free module with a single generator v\ and relations 



EiV\ = 0, KiV\ = q 



Remark. Here and below we work over the field F = C(q)({q a , a G C}). In this 
setting, g^' A ) is a function ()* — > F. 
We have the decomposition 



M x = Ma [A — /x] 
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of M\ into direct sum of weight subspaces M\ [A — ji] , where we say that a vector 
v has weight /j, G f)* if 

= g^^r; 

The restricted dual module M^ is a W + -module with a lowest weight vector v*_ A 
such that (v* A , u a) = 1. By definition we have 

(V» = ^ G M A ,v* G M A *. 

Introduce a symmetric form Fa on Ma defined by 

F(giV X ,g2V\) = (w(5l)v* A ,02VA)- 

The weight subspaces are pairwise orthogonal with respect to this form. The re- 
striction of F to weight subspaces M\-^ is proportional to the quantum Shapovalov 
form J 7 , introduced in [CK]: 

for some constants C M . 

Fix a basis G U~[/j]. Let F M = (-F M )ij, i, j = 1,2,..., dim Ma [A — /j], denote 
the matrix of the restriction of form F to Ma [A — fj] with respect to the basis gtfvx 
of Ma [A — fj] . A variation of the quantum determinant formula [CK] asserts that 



detF M (A)=c n Yi{i-q- 2 ^ x+p)+n{a ^y ar{fi na) 



where Par is the generalized Kostant partition function, and C G C(q) is a constant, 
depending on the choice of basis gf. 

This determinant is a linear combination of terms q~ 2 ^ ,x \ where /it's belong to 
a finite subset L C Q, with some coefficients from C(q). Motivated by this fact, we 
introduce 

Definition. Expressions of the form X^eL a^q -2 ^^, G C(q), will be called 
q-polynomials with support L and coefficients a M . 

Verma modules are reducible when form F is degenerate, i.e. detF M (A) = for 
some ii. This happens when A satisfies one of the Kac-Kazhdan equations: 

Tl 

(2-1) (a,X + p) = -(a,a), n = l,2,... 

For A generic from Kac-Kazhdan hyperplanes, Ma contains a unique submodule 
M A , isomorphic to M\- na . 

Intertwining operators. Let U be an irreducible finite-dimensional U q g-module 
with non-trivial zero weight subspace U[0]. For u G U let $ A : M A — > Ma <E> U be an 
intertwining operator such that i>a -^Da®« + higher order terms, where "higher 
order terms" mean terms of the form V\-^ ®u^,ijl > 0. 

If Ma is irreducible, then $ A exists and is unique for any u G U[0]. Indeed, 
we have a unique W + -intertwiner O : M A — > E7, such that Ot> * A = it. Since 
Hom(M x , U) = M A *<g>?7 = M A <g>£7, it corresponds to a singular (i.e. W + -invariant) 
vector <fi G Ma <8> We now construct $ A by putting $ a ^a = and extending $ A 
to the whole Ma by the intertwining property. 

For our purposes we need an explicit form for that singular vector. 
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Proposition 2.1. For any (homogeneous) basis {gf} o/W 

(2-2) = e [Eft" 1 ),- (%r«A ® WW 



it 

?J 



is a singular vector in M\ <g> £/. 

Note that since U has a highest weight, the summation is over the finite set of 
/Li's such that ^ 0. 

Proof. We check that the corresponding element $ G Hom(M£, U) defined as the 
composition — > <g> Ma <8> ?7 — > U is a W + -intertwiner. We have: 



* M<7>Ia)) = E E Ms>* a,<7>aM<7> = 

E (*7% WKa,J?«aW)« = E ( E Wni) = 

i,j j \ i J 

J2$3nU(gj)u = U}(gn)u 
3 

□ 

Recall that in classical case matrix elements of the inverse matrix F~ were 
rational functions of A with at most simple poles in the Kac-Kazhdan hyperplanes 
given by (2-1). (see [ES]). A similar argument, also involving Jantzen filtration, 
proves that the same is true in the quantum case. Therefore, expression (2-2) for the 
singular vector is a ratio of two g-polynomials, with at most simple singularities 
in a finite collection of Kac-Kazhdan hyperplanes. 

If we multiply the (/-rational expression (2-2) by the least common denominator 
x(A), we will get a well-defined for all A's formula for a singular vector G M\ ® U. 
We are now going to show that in fact the least common denominator may only 
contain factors 

X «(A) = l_ ? -2<a,A+p>+n< a ,a> 5 

corresponding to n, a such that [/[net] ^ 0. Indeed, suppose that x(\) contained a 
factor x"(A), but [/[net] = 0. 

Consider A generic from the hyperplane (a, A + p) = ^(a,a). Then M\ contains 
a unique maximal submodule M\ = M\_ na , generated by singular vector v\— nOL . 
Since the first term x(ty v \ ® w in the expression for turns into zero in our 
hyperplane, the singular vector must have form 

4> = V\- nct <S> u + higher order terms 

The intertwining property implies that u G (7[na], and by assumption u = 0. 
Therefore, <p is zero f° r A generic from the hyperplane, and by Bezout theorem is 
divisible by Xn- This shows that x(\) was not the least common denominator - 
contradiction. 

Denote k a = ma,x{n\U[na] ^ 0}, 

Lg = \y\v = E] rn a a, < m a < k a for all a G R + } 

aER+ 

We conclude this section with the following 
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Proposition 2.2. IfU is an irreducible finite- dimensionalU q g-module with highest 
weight 6, the singular vector G M\ <g> U, given by (2-2), can be represented as 

, Ez Si(\)giv\ g ui 

(1 - g -2<a.A+P>+m(a,a>) 

where gi G U~,ui G ET, and q-polynomials Si(X) have supports, contained in Lg. 

Proof. We already proved that the least common denominator for the expression 
for (f) may only contain factors %"(A), n = 1, . . . , fc a . 

The statement about the support of the polynomials Si(X) follows from the fact 
that the support of the numerator must lie within the convex hull of the support 
of the denominator, which in this case is exactly Lg. □ 

3. Matrix Trace, the \l/-function and its properties. 

We now fix an irreducible U q Q-modu\e U with highest weight 6 and non-trivial 
zero weight subspace. We use the notation 

k a = ma,x{n\U[na] ^ 0}, = k a • a G Q+ 
XZ(X) = 1 - q- 2{a ' X+f>}+n{a ' a) , X(A)= J] n X £(A) 

aeR+ n=l 

Lg = {p G Q+|/i = m a a, < m a < k a for all a G R + } 

aER+ 

As in [ES], define a new intertwining operator 

^ = x(A)^:M A -M A ®tf 

From Proposition 2.2 it follows that is well-defined even when A belongs to 
Kac-Kazhdan hyperplanes, where did not always exist. 
Introduce a End (U[0]) -valued function *(A, x), A, x G f)*, by 

*(A,^ = Tr| MA (^e-) 

Proposition 3.1. The function ty(X,x) defined above has the form 

®{X,x) = e < A ' x > q~ 2{lI ' X+p) P»(x), 

where P^(x) G End(U[0]) and P&(x) is invertible for generic x. 
If (a, A + p) = §(ck, a) for some a G R + , n = 1, 2, . . . , k a then 

(3-1) V(X,x) = ^(X-na,x)B nce (X) 
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for some (possibly infinite) sum 

B na {\) = T 2<M ' A> ^ e End(U[0]) 

(In fact, the matrix elements of B na (X)'s are ratios of q-polynomials.) 

Remark. If we take U to be a trivial module, then the \I/-function becomes the 
usual character of the Verma module. Therefore, we can regard the ^-function as 
a generalized (matrix- valued) character of the Verma module M\. 

Proof of Proposition 3.1. 

Tr\ Mx {®le*)=e^ £ e -^>B M (A)«, 

where "partial traces" £> M (A) € End(U[0]), corresponding to weight subspaces 
M\ [A — p] , are defined by 

£ M (A)w = Tr (Proj MA[A _ M] o Proj MA[A _ M] ) 

By Proposition 2.2, <B M (A) are End (U[0]) -valued (/-polynomials with support Lq. If 
we let x — > oo in such a way that — > +00, ct e -R+ (this just means that we 
are keeping only the highest weight terms of the series) , we will get asymptotically 

Therefore Pq{x) ~ 1, and P®(x) is invertible for generic x. 
We now prove the second property. 

If A is generic from hyperplane (a, A + p) = ^(a,a), then M\ is reducible and 
contains a unique submodule M\ generated by singular vector v\- na . It is clear 
that for such A there are no order p terms in 3>™v\ unless p > na. In other words, 
$ A maps M\ into M\ <g> U, and 

— v \— na ® u + higher order terms 

$\V\- na = v\- na <g) u" + higher order terms 
Clearly, both u' and u" depend linearly on u. More precisely, 
(3-3) u" = B na {\)u, 

Therefore we have 

X {\-nam\ Mi =n_ nct 

Taking the traces of the operators from the last equation and using (3-3), we get 

X(A - na)*(A, x) = *(A - na, x)B na (\) 

Since %(A — na) is invertible in the "Laurent series" completion of C[P], we can 
introduce 

B (\) - B ^( A ) - V a 2(tl ' x) B^ 
y( A — na) ^— ' 

and (3-1) follows. Proposition 3.1 is proved. □ 

We prove that property (3-1) of the \l/-function determines it uniquely up to 
multiplication by a factor, depending only on x. 
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Proposition 3.2. Suppose we have an End(U[0]) -valued function 



*'(A, x) = e^ Q~ 2ifl,X) Q^) 

where L C Q+ and Q^(x) G End(U[0]), satisfying condition (3-1). Then the L 
contains at least one weight \i > ©. 

Proof (cf. [ES]). Let us rewrite the condition (3-1). We have: 

For every /i G Q consider the set 

Z a (n) = jU + = {z/GQ|z/ = / u + ma for some m G Z} 
Comparing coefficients for q^> x \ we get for n = 1, 2, . . . , fc Q : 

This is a system of linear equations on unknown functions P u (x). Note that the 
summation over [3 is finite, since P v -p = if v — f3 ^ Q_|_. The matrix of this 
system has a block-upper-triangular form, blocks corresponding to subsets Z a (/j,) 
for different \i. The determinant of this matrix is an entire function of x, and 
asymptotically when (a, x) — > +oo, a G -R+ 

lhs ~ en{a,x) q~ n{ " ,a) 'Pu(x) 

Therefore, asymptotically the determinant of the matrix of this system is equal to a 
Vandermonde-type determinant, which is nonzero. It follows that the determinant 
of the system of equations is nonzero for generic x. 

Suppose \i is such that not all P^x) is not identically zero (i.e. we have a 
nontrivial solution of the system of equations). Then for such fx we need to have 
Card (Z Q (//)) > k a for all a. 

The rest is based on the following combinatorial Lemma. 

Lemma 3.3. If a subset L C Q+ is such that for any a G -R, \i G L we have 

(3-4) CardfX^HL) > k a 

then any maximal weight fxo G L satisfies 

Ho > Q = ^ k a ■ a 

OL 
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Proof of Lemma. We start by proving that (// ,Ai) > (0,Ai). To any weight \x 
assign a string of numbers 

{(//,Ai),...,(//,A r )} 

where A l5 A 2 , . . . , A r are the fundamental weights. Take the set of all roots f3 G R + 
such that (/?, Ai) > and arrange them in the increasing lexicographical order: 

01 ~< (3 2 ~< < 0N 

Our goal is to construct a "descending" sequence of weights /ii, . . . , jiN G L such 
that 

(3-5) Ai) - (//,-, Ai) > /c ft Ai) 

Then we conclude that 

N 

(3-6) {fio, Ai) > (/Mr, Ai) + ^ fcp. Ai) > + ^ fc a (a, Ai) = (0, Ai) 

i=i aefl+ 
We now explain how to construct jttj+i given //j. 

We will build an auxiliary sequence fip i = 0,1,2,...,^ such that three 
conditions are satisfied: 

r 

(*) jU* + n/9j + ^ niai G L, n, > 0, implies all ni = 

(**) + nj3j G L for at most % different n G N, 

(***) (^ +1 ,A 1 )<(^,A 1 )-(/3 J ,A 1 ) 

We take $ = fij. Suppose we have already constructed /i*-, i < k^, satisfying 
(*) and (**). Conditions (3-4) and (**) guarantee that we can find in L a weight 
of the form — n(3j, n G N. Take the minimal positive integer m such that the set 

{/j G L\ (fj,, Ai) = - mPj, Ai), // y fi*j - m(3j} 

is nonempty, and let \x^ x be its maximal element in the sense of lexicographical 
order. We claim that satisfies (*), (**) and (***). Indeed, (***) is clear from 
the construction. To prove (*) we note that if we have 

r 

H = + n(3j + ^2 n i a i e L 

1=2 

for n, ni > such that not all of them are equal to zero, then either n > 1 and m was 
not minimal, or n = and (*) was not satisfied for - contradiction. Finally, (**) 
is obvious when — = mf3j] if it is not the case, i.e. /ij — = mPj a h 
then there is no n G N at all such that \i l ^ x + n/3j G L. 

Having thus constructed all fip i = 1,2,..., kp j , we can put fij+i = [i j 3 ■ From 
construction we see that (3-5) is satisfied, and by (3-6) we conclude that indeed 
(fio, Ai) > (0, Ai). We can prove similarly that 

(/io,Ai) > (0,Ai), z = l,2,...,r 

Therefore, /j, > 0. The Lemma is proven. □ 
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Corollary 3.4. The \& -function, satisfying (3-1), is unique up to a factor, depend- 
ing on x. 

Proof. If we have another function ^'(A, x) with highest coefficient Pq{x), satisfying 
(3-1), then the function 

0(A, x) = *'(A, x) - P'®{x) (Pe(x))- 1 *(A, x) 

will still satisfy (3-1), but its support will only contain weights \x < 6. By Propo- 
sition 3.2, 0(A, x) = 0, and the statement follows. □ 

Proposition 3.5. Any function (f>(X,x), satisfying (3-1), which has the form 

<f>{\x) = e^Q(X,x), 

where Q(A, x) is an End (U[0]) -valued q-polynomial with depending on x coefficients, 
can be represented as 

4>(X, x) = c(A, x)ty(X, x) 

for some End(U[0]) -valued q-polynomial c(X,x) with depending on x coefficients, 
satisfying 

cl^X-—,x) = c[X+ —,x) 

The proof uses induction on support of c(A) and is very similar to the one in 
[ES]. However, we will not need this fact for our main theorem, and we leave the 
details to the reader. 

4. Existence of difference operators. 

Introduce a family of difference operators T\, corresponding to weights A e P, 
acting on ()*, by 

T A (x) =x + Alogg 2 , x E fj* 
They naturally act on functions on ()*; for example, for f{x) = e^ a,x ^ we have 

(T A f)(x) = e (^+ Ai °Sq 2 ) = ? 2(a,A) e <«,x> = q 2{a,K)j^ 

Proposition 4.1. Any function (f)(X,x), satisfying (3-1), which has the form 

(P = e^ x >P(X,x) 
for some q-polynomial P(X, x), can be represented as 

<p(X,x) = D^(X,x) 

for a unique difference operator D with depending on x coefficients. 

Proof. Uniqueness is obvious, since otherwise the ^-function would be annihilated 
by a nontrivail difference operator for all A, which is impossible. (See, for example, 
[EK]). 
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To prove existence of D we use induction on the support of P(X,x). Consider 
the family of all finite subsets L C Q such that for /i, v G Q 

Such subsets are all finite unions of cones of the form C M = \y G Q|z/ < /u}. 

We prove that if our statement is true for all (/-polynomials P(X, x) whose support 
is strictly contained in L, then it is also true for (/-polynomials with support L. 

If L does not contain weights \x > 6 then by Lemma we have P(A, x) = 0, and 
the operator D = 0. 

Suppose there is a z/ G L such that i/ > ©. Consider the set of all such z/s and 
let // be a maximal element from this set. 
Consider the function 

(f>'(\, x) = 0(A, x) - P M (x)T M _ e * (A, x) 

It satisfies (3-1) and has support strictly contained in L. By induction hypothesis 
we can represent 4>' as 

0'(A,x) = D'*(\,x) 

The operator 

D = P^{x)T^q + D' 
satisfies the required properties. □ 
We now prove a simple technical 
Lemma 4.2. Let a q-polynomial c(A) be represented as 

c(A)= £ cv(A) 

ttEP/Q 

where c^{\) are q-polynomials with support in the coset tt + Q. 
Suppose for some n, a we have that 

(4-D c ( A+ _)=c(A--) 

whenever (a, A) = 0. Then (4-1) is also satisfied for each c^{\). 

Proof. Property (4-1) is equivalent to divisibility by (/ 2 < a > A > — 1 of the (/-polynomial 

c(A) = c(A + T )-c(A- T ) 

On the other hand, one can see that 

c(A)= 5 ^ A )' 

7T€P/Q 

where 

Ctt(A) = c n \\ + — j - c n (A - — J 

Clearly, c(A) is divisible by (/ 2 ( a > A > — 1 if and only if all c,r(A)'s are divisible by 
q,2(a,A) _ ^ an( j Lemma follows. □ 
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Theorem 4.3. For any q-polynomial c(A) such that for any a G R + 

. « /, na\ / na\ 

(4-2) c (a+— )=c(A-— J, ra = l,2,...,fc a 

whenever (a, A) = 0, t/iere exists a difference operator D c with coefficients in 
End(U[0]) suc/i <fcat 

L> c *(A,x) = *(A,x)c(A + p) 

The correspondence c(A) — > _D C is a homomorphism of rings. 

Remark. We put c(A+p) on the right since in that form it admits generalization 
to the matrix case (see Theorem 4.4). Of course, for scalar (/-polynomial c(A) we 
could write it in a more traditional form D C \I/(A, x) = c(A + p)\I/(A, x). 

Proof. By Lemma 4.2, it suffices to prove the theorem for c n (X) of the form 

C (A) = q 2 ^ c^q- 2 ^ 

for some fio G P. Consider the function 

<t>{\x) = T_^(X,x)c(X + p). 
It satisfies (3-1), and it has the form 

<t>(X,x) = e^ Q-^QM, 

By Proposition 4.1 it can be represented as 

<f>{\x) = DV(\,x) 
for some difference operator D. Put D c = T^ D. Then we have 

D c tf(A, x) = T Mo Dtf(A, x) = T Mo 0(A, x) = 
T M0 T_ M0 *(A, x)c(A + p) = *(A, x)c(A + p), 

We now prove the homomorphism property. Suppose we have two polynomials 
c(A) and c'(A). It is easily checked that operator D cc > — D C D C > annihilates the *f>- 
function for any A; therefore it has to be identically zero. Hence our correspondence 
is a homomorphism of rings. □ 

We have big supply of (scalar) (/-polynomials, satisfying (4-2), arising from the 
algebra of Weyl group invariant (/-polynomials, which is freely generated by the 
Casimir elements c\,...,c r . This give us n algebraically independent difference 
operators Di, . . . , D r . However, there exist other (/-polynomials, with this property. 
For instance, 

k 

c (a)= n n (q 2{ ^ +n ^-i) 

a6i?+ n= — k a 

also satisfies (4-2). It gives rise to a difference operator, commuting with all those 
generated by the Casimir elements, but not lying in the ring generated by them. 
This procedure gives examples of what we called algebraically integrable commu- 
tative rings of difference operators. 

Theorem 4.3 can be slightly generalized to the matrix case. 
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Theorem 4.4. 

For any End (U[0]) -valued q-polynomial C(X) such that 

C ( X+! Y) Sna(A) = B ™W C ( A - "y) ' aeR+,n=l,...,k a 

whenever (a, A) = 0, there exists a unique difference operator Dq with coefficients 
in End(U[0]), such that 

D c $(X,x) = ^(X,x)C(X + p) 

The correspondence £ : C(A) h- > -De a homomorphism of rings. 

Proof. The argument used in proof of Theorem 4.3 in the obvious way extends to 
the matrix case. □ 

Remark. Operators Di, . . . , D r act on \l/-fuction as scalars, and therefore com- 
mute with all operators Dq constructed as above. In fact, one can show that the 
centralizer of the subring generated by operators Di, . . . , D r in ID£ (?/[()]) coincides 
with the image of £. We do not include the proof of this statement here. 

In the next section we explain how our construction is related to Macdonald 
theory. 



5. Root system A n and Macdonald theory. 

Consider a special case of our construction for q = sIn, U = Uk = S^ kN ^V, 
where V is the fundamental representation. It is well-known that the zero weight 
subspace U[0] is one-dimensional, and we can regard ^(A, x) as a scalar-valued 
function. Note also that in this case k a = k for all a G R + , and = k J2 a eR+ a = 
2kp. We have di = 1, and therefore qi = q for alH = 1, . . . , r. 

We first prove an important property of partial traces 

B M (A) = Tr (Proj MA[A _ M] o Proj MA[A _ M] ) , 

introduced in Section 3. We will use the notation 

q n_ q -n q 2n -l l-q- 2n 



q — q q — 1 1 — q 

Proposition 5.1. Given any a G -R + , n = 1, . . . , k, we have for all n e Q 

(5-1) B M (A) = q~^ a ^ ■ S M _ na (A - na) 

whenever (a, A + p) = n. 

Corollary 5.2. For a G -R + , n = 1, . . . , k, the function \&(A, x) satisfies 
* (A + f,x )=,-<««» (A -™ x) 

whenever (a, A) = n. 

Proof of Proposition 5.1. For given n, a it is sufficient to prove a special case of 
(5-1), corresponding to p = na : 

(5-2) B nce (X) = q-^ a ^ ■ B (X - na) = (T <na,e> x(A - na) 
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Indeed, let A be such that (a, A + p) = n. Then the image of is contained in 
M\_ na ® U. From (5-2) we see that 

and the more general formula (5-1) follows. We now use induction on the height of 
root a to prove formula (5-2), and thus the Proposition 5.1. 

Base of induction. Consider the case when a = oti is a simple root. Let A 
be generic in the hyperplane (aij, A + p) = n. Then M\ contains a unique nonzero 
proper submodule M\, generated by the singular vector v\- ncti = F™v\. One can 
check that 



L J + * \m=n+l / Wa,m=l 



^A-na ® + t.O.W 



where t.o.w. denotes "terms of other weights" (in the first component). 
It follows that 

*wa) = {q ~ l -f ( n (i - q - 2{n - m) )\ [ n n x&w ) w 

L J + ' \m=n+l / \^ aj m=l / 

[fc+n 



It is known that F^E™ acts as multiplication by |^^ | in ?7[0] . Since 



II 



L J+ \m=n+l / L J m=l 

for A from the hyperplane (ojj, A + p) = n we have 

^(A) = q -**» n (i - ? 2 ^) ( n n 



m=l \ (3^a.i m=l 



q~ 2kn II XS? (* - nil^ A - ^) ) = Q'^^xiX - not) 

\m—l J V 7^a, m=l 



Induction step. Suppose (5-1) is true for all roots (3 such that height (3 < 
height a. We are going to prove that (5-2) is true also for a. 

We first show that B na {\) is divisible by factors x^(A — na) for all f3 ^ a, m = 
1, . . . , k. It suffices to prove that B na (X) vanishes whenever (/?, A — na + p) = m. 

Consider two cases: 

(1) If 8 a (J3) G R + , put 7 = s a ((3). Then (7, A + p) = (0, A - na + p) = m. 
If A is generic from hyperplane (7, A + p) = m, then the image of 
is contained in M\- mi . Since not — 7717 ^ Q+, there will be no terms con- 
tributing to B na {\) = Tr \ Mx[x _ nce] ^ and ^na(A) = generically (and, 
therefore, identically) in the hyperplane (7, A + p) = m. 
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(2) If s a ((3) £ R + : put 7 = —s a {(3) E R + . Then a = f3 + 7, and we can assume 
the induction hypothesis true for (3 and 7. 
We have: 

(7, A + p) = -(/?, A - na + p) = —m 

(/?, A + 7777 + p) = (P,X — na + p) + n((3, a) + m(7, (3) = m + n — m = n 
Therefore, 

B na (X) = q^^B na+mi {\ + mi) = q^- n ^B (n+mh (X + m 7 - n(3) 
Also, 

(a, A + rwy — n(3 + p) = n + m — n = m 
so the image of is contained in M A _( m+n ) /3 , and 

B {n+m)l (X + mi - n/3) = 

because ma — (m + n)(3 Q + and there are no terms contributing to 

B {n+m)l (X + mi - n(3). 

We have proved that B na (X) vanishes in the required hyperplanes, and is there- 
fore divisible by all the required factors. Thus, in the hyperplane (a, X + p) = n we 
get 

k 

B na (X) = C(X) J] \[ X i{X-na) 

for some (/-polynomial C(A). It is easy to see by comparing highest terms that 
C(A) is constant in the hyperplane (a, X + p) = n. To compute this constant, take A 
generic such that (a, X + p) = (/?, X + p) = n. Then automatically (7, X + p—n(3) = n. 
We have: 

*(A, x) = q- {nf3 ' e ^(X - n/3, x) = g-<^e) ? -<n 7 ,e)^ (A _ n(3 _ n ^ x ) 

But A — n(3 — ni = X — na does not lie on any Kac-Kazhdan hyperplanes, therefore 
\& (A — na, x) 7^ 0, and 

Proposition 5.1 is now proved. □ 

Let A be a dominant integral weight, and V\ is the irreducible W g g-module with 
highest weight A. 

Proposition 5.3 (Generalized Weyl formula). The operator $™ : M\ — > M\® 
U descends to a homomorphism V\ — > V\ <g> U. The function 

p x (x) = Tr\ Vx ($ u x e x ) 
is expressed in terms of the functions ^(A, x) by 

(5-3) q {e > x) p x (x) = ^(_i)'W g (eyA) $(w p A)j) 

wew 
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Proof. The operator defines an operator M\ — > V\ ® U. This operator has to 
factor through V\ because it lands in a finite dimensional representation. Thus, $^ 
in fact defines an operator V\ — > V\ <S> U. 

Recall that for A G P++ we have a resolution 

0^V x ^M° x ^Ml^M 2 x ^..., 



where 



= M\, 

For matrix traces we have as for usual characters 

p x (x) = Tr\ Vx me*) = £(-l)*Tr|^ (^e a 



When A is generic from hyperplane (ctj, A + p) = ^(a^, ctij), n > k, then Ma 
contains a submodule Mj[ = M\- nai , generated by a singular vector v\- nai = 
F™v\. Since Ma [A — ncti] is one-dimensional, we can write 



. = Vx-n^ ® + • • • 



Then we will have ^aIm 1 = ^A-n ai - To compute it" we use the formula 



nvx = E (9 ~ N+ f ( n xr(A) ) ( n Ilxfw ] F>A®i?>+ 



m=0 



v 2=m+l 



/3^ a< Z=l 



where we retained only terms, which will contribute to the expression for u" . Ap- 
plying A(F™) to the RHS, and collecting terms involving v\- nai , we deduce that 



m=0 



n 
m 



(g-'-g)' 

[m]+! 



II ^"(A) n nxf(A) br«r- m sr« = 



E(-i) m g m_fc 



m=0 



n 
m 



l=m+l 

(g _1 — g) fc [n — m — 1]_! 
[to] + ! [n-fc-l]_! 



v /3^a t Z=l 
fc 

n n^( A ) i g 2m(n_m) ^r^> 



The operator F^E™ acts in £/[0] as multiplication by 
Also, we need to use the identity 



[n]! 



[n-m-1]! [fc + m]! , + 



m=0 



[m]![n-m]! [m] ! [n - - 1] ! [fc-m]! 



= (-1)' 



[n]\ ■ 



It can be interpreted as equality of two polynomials in z = q 2n of degree k. To 
prove this identity it suffices to check it for n = 0, . . . , k, when there is only one 
nonzero term in the LHS of the equation. 
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After easy transformations, we conclude that 

k 



u" = q~ 2kn 



n a - <r 2(n+m) ) n n^(A)i«= 

\m=l / \[3=£ ai m=l 



q~ 2kn ( n^(A-naO] (il II *-( A " J M = 9~ n<a *' e> x(A - naju 

\m—l / V -f^cti m=l / 

It follows that the restriction of $^ to the submodule Mj[ = M\_ na . coincides with 
q~ n ( ai ' e '®\- nai for A generic from hyperplane (a^A + p) = ^(cn^aij). Therefore, 
it is true for all A from that hyperplane. Even more generally, for A G P++ we have 
4 <e ' A> *AVp A = q {e ' wPX} K PX - Thus 

q {6 > X) Tr\ Mi ($ u x e*)= £ 4 <e ' A> Tr| M ^ A (^e*) = £ g< e ^ PA >*KA,x) 

l(w)=i l(w)=i 

Therefore 

□ 

Introduce a normalized matrix trace by 

g< G ' A -^(A -p, x) 



Condition (3-1) can be rewritten for the function ip(\,x) as 
(5-4) ,( A + !£) 

for all a G -R + , n = 1, . . . , k and for all A such that (a, A) = 0. 

Corollary (Macdonald polynomials). Combining (5-3) with results in [EK], 
we see that (up to a factor) Macdonald polynomials are equal to 

(5-5) p x (x)= ^2(-l) l ^(w(X+(k + l)p),x) 

wew 

Proposition 5.4 (Macdonald operators). The function x) is the common 
eigenfunction for Macdonald operators, corresponding to t = q k+1 : 

I k+l e (oi,x)/2 _ -fc-1 -<a,x)/2 \ 

-^=5> II e (^)/2 _ e -(^)/2 T Ai , i=l,...,r 

\(a,A i ) = l / 

TTie corresponding eigenvalues are W -invariant q- polynomials q(A) = J2 weW q( X ' wAi ) 

Miip{\,x) = Ci{\)ip{\,x) 

Proof. We already know that the \l/-function is the common eigenfunction for a 
family of commuting difference operators, corresponding to c(A), satisfying (4-2). It 
is shown in [EK] by using central elements in U q Q, that for the \l/-function normalized 
as above, the operators, corresponding to elementary symmetric functions q(A), are 
exactly Macdonald operators. □ 

Now we study relations between unctions for different values of k. We use 
notation ipk(\,x) for the •i/'-function constructed from representation Uk- 
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Theorem 5.5 (Shift operators). There exist difference operators G k such that 

ip k+1 (X,x) = G k ip k (X,x) 

These operators are W invariant, and there action in the basis of Macdonald poly- 
nomials is given by the formula 

G k Pk,i( x ) = Pk+i,i- P (x)J -pe P+; G k p k ,i = 0, I - p £ P + . 

Remark. Shift operators in the q-deformed case were introduced by Cherednik, 
[Chi]. 

Proof. The argument from the proof of Proposition 4.1 can be used to prove that 
any function, satisfying (5-4), can be represented as 

0(A,x) = Dip k (\,x) 

for some difference operator D. Applying this to the function ^ +1 (A, x), we prove 
the existence of an operator G k such that 

ip k+1 (X,x) = G k ip k (X,x) 

From the generalized Weyl formula we get 

G kPk , x = (-l)' (w) GWfc(«>(A + (k + l)p), x) = 

wew 

E (-^) l(w) Mw(X + (k + l)p),x) = Pk , x - P (x) 
wew 

If I — p is not dominant then the right hand side of the last formula is zero by (5-4). 

We see that G k maps Macdonald polynomials to VF-invariant functions. This 
implies that G k is itself VF-invariant. □ 

Remark. We saw that shift operators relate eigenfunctions for Macdonald op- 
erators, corresponding to different (integral) values of parameter k. One can write 
it in the form 

(5-6) G k M\ k) = M\ k+1) G k 

One can check that G k analytically depends on k, therefore one can extend equality 
(5-6) to the case of arbitrary k. This implies the existence of shift operators in the 
general case, which is proven in [Chi] using representation theory of double affine 
Hecke algebras. 
Denote 

s k (x)= J] fj (y e <^>/ 2 - <rV<^>/ 2 ) 

aei?+ i=-k 
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Theorem 5.6 (Duality). The function <pk(\ x ) = ^(^OV^A, x) is symmetric 
with respect to transformation q( a > x } <-> e^ a ' x \ 

Proof. The idea of the proof is the same as in [VSC]. First we prove that (pk 
as a function of x satisfies condition (4-2), and then duality will follow from the 
uniqueness property. 

We already know that ipk is the eigenfunction for the Macdonald operators. 
Therefore, (fk{\ x) is the eigenfunction for the operators, obtained from Macdonald 
operators by conjugation by Sk(x). Such an operator, corresponding to a minuscule 
weight A, is also VF-invariant and has the form 

M A = /aH t »(A). 
wew 

where 

q -k e {p,x)/2 _ q k e -{p,x)/2 

AO) = || e {/3,x)/2 _ e -{(),x)/2 

(/3,A> = 1 

Fix a root a and denote W a = {w G W\ (a, wA) = 1}. Let x be such that (a, x) = 
0. The </>-function does not have singularities along that hyperplane. Collecting all 
the singular terms in the equation 

M^k{\ x) = c A (\)tp k (\, x), 
which occur forweW such that (a, wA) = ±1, we obtain 

(5-7) Y f^( wx ) T w{A)^k(\ x) + Y fA(s a u)x)T (SaW){A) ip k (\, x) = 

wew a wew a 

When x belongs to hyperplane (a, x) = we have 

q -k-l e (P,x)/2 _ q k+l e -{(3,x)/2\ ( q -k e {/3,x)/2 _ q k e -{f3,x)/2 

e (/3,x>/2 _ e -(f3,x)/2 \=-(s a w)\ || (/3,x)/2 Z e -(/3,x)/2 

/ \</3,A) = l 

therefore (5-7) simplifies to 




^ /a(wx) (ipk (A, x w + ^ - y? fc ^A, 

-i 

where 



x w -^)) = 



= x + w(A)" 1 " = x + w(A) — — 



also belongs to hyperplane (a,x w ) = 0. It follows that 

<Pk (A, x + ^ - p fc (\ x - = 

identically when (a,x) = 0, i.e. that condition (5-4) is satisfied for n = 1. Taking 
n-th power of operators .Ma, and repeating the same argument, we can prove that 
it is satisfied for n = 1, . . . , k. 
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From the obvious modification of the theorem about uniqueness of the \I/-function, 
we conclude that the function (pk(X,x) transforms into itself when we interchange 



Remark. This duality result is closely related to the symmetry of the difference 
Fourier pairing defined recently by Cherednik [Ch2]. 

Appendix. Example: g = sl 2 . 
In this case we have only one root a, and we can make identifications 

e {X,x) ^ e \x^ e («,x> ^ e 2x ? Q {a,X) ^ g X 

Operators T, T~ l act on function f(x) by 

(Tf)(x) = f(x+1), (T- 1 f)(x) = f(x-l) 

Case k = 0. This is the simplest example, corresponding to the case U - trivial 
representation. In that case 

$lv x = ®\v\ = v x ®u 

Xx (\+l)x 



1 - e~ 2x 



— e 



e Xx 



— e 

Formula (5-5) is the usual Weyl formula 

e (A+l)x _ e -(A+l)x 

char V\ = 

e x — e x 

The operator, corresponding to VF-invariant polynomial ci(A) = q x + q~ x , is the 
Macdonald operator 

_ qe x - q- 1 e~ x T qe~ x - q~ x e x T _ x 
and we have 

.MiVu(A, x) = ci(A)^i(A, x) 

Condition (4-2) gives no restriction on c(A), and we can set co(A) = q x . The 
corresponding operator Mo, whose existence is predicted by Theorem 4.1, is equal 
to 

M = 9eX ~ 9 " le " X T, 
e x — e x 

and Mi = M + Mq 1 . 
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Case k = 1. Let U now be the 3-dimensional representation. We have: 

1 — q~ IX 

= (! - 9~ 2A )^A ®u-(q- q~ 1 )Fv x <g> £w 

„Ai / „— 2x 

*(^> = T^(i-<r 2A -(<; 2 -<r 2 ) r 



g 2 e 2x 



g Ax f 1 — q 2 e 2x 



ipi{\x) = 



q~ 2X 

1 — e 2x \1 — q 2 e 2x 



e x _ e x yg e x _ q i e x ^ i e x _ ^ e x 



It is easy to check that indeed 

*(0, x) = Q" 2 ^(-2, x), x) = Vi(-1, x) 

Function ^(A, x) is related to ipo(X,x) by 

V>i(A,x) = G ^o(A, x) 
where the shift operator Go is equal to 

G = ^— (T-T- 1 ) 

The Macdonald operator .Mi is equal to 

_ q 2 e x - q- 2 e~ x T q 2 e~ x - q~ 2 e x 

.MiV>i(A,x) = ci(A)V'i(A,x) 

The operator .Mo, corresponding to the eigenvalue c (A) = g 3A — [3] g (/ A , is equal to 

q 4 e x — q~ 4 e~ x q 3 e x — q~ 3 e~ x 3 ge x — (/ _1 e _x 
Alo = ~i — — — T — |3L^— : ^— T 

qe x — q e x e x — e x q e x — qe x 

Classical limit q — >• 1. Set e = log(q). We have the expansion 

V>i(A,x) = 2e-Vi 0) (A,x)+O(e 2 ) 

where 

, % „Ax / „x I „ — x 

V 0> (A,x) = - 6 (a-V^ 
( e x _ e -x) z y e x — e x 

is the •i/'-function for the classical case (cf. [ES]). The difference operators become 

Mi = 2 + e 2 (T? 2 + 4) +0(e 3 ) 

M = -2 + e 2 (3T? 2 + 8) + 4e 3 (P 3 + 1) + 0(e 4 ), 

where commuting differential operators P 2 ,^3 are equal 

<9 2 e x + e _x 9 
^2 = -tt^ + 4- 



<9x 2 e x — e x dx 
d 3 fe x + e~ x \d 2 ( 12 \ d /e 3x - 3e x - 3e" x + e" 3i 



<9x 3 \e x — e x ) dx 2 \ (e x — e x ) 2 J dx \ (e x — e x ) 3 
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